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Abstract 


GLACE is a model intercomparison study focusing on a typically ne- 
glected yet critical element of numerical weather and climate modeling: land- 
atmosphere coupling strength, or the degree to which anomalies in land sur- 
face state (e.g., soil moisture) can affect rainfall generation and other at- 
mospheric processes. The twelve AGCM groups participating in GLACE 
performed a series of simple numerical experiments that allow the objec- 
tive quantification of this element. The derived coupling strengths vary 
widely. Some similarity, however, is found in the spatial patterns gener- 
ated by the models, enough similarity to pinpoint multi-model “hot spots” 
of land-atmosphere coupling. For boreal summer, such hot spots for precipi- 
tation and temperature are found over large regions of Africa, central North 
America and India; a hot spot for temperature is also found over eastern 
China. The design of the GLACE simulations are described in full detail so 
that any interested modeling group can repeat them easily and thereby place 
their model’s coupling strength within the broad range of those documented 
here. 
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1 Introduction 

1.1 Motivation 

Precipitation has a clear impact on soil moisture: large rain events tend to 
wet the soil. To what extent, though, do land surface moisture and tem- 
perature states affect in turn the evolution of weather and the generation of 
precipitation? How does a human-induced change in land cover affect local 
and remote weather, if at all? Such questions lie at the heart of much recent 
climatological research. This research is largely performed with numerical 
models of weather and climate (atmospheric general circulation models, or 
AGCMs, and mesoscale models), mostly because direct observations of the 
impact of land surface anomalies on atmospheric behavior are difficult, if not 
impossible, to obtain at regional to continental scales. Also, AGCMs and 
mesoscale models have the advantage of being amenable to sensitivity stud- 
ies - their process parameterizations can be manipulated easily in controlled 
experiments. 

The list of published AGCM studies that address questions of land- 
climate interaction is extensive (e.g., Charney et al., 1977; Shukla and Mintz, 
1982; Henderson-Sellers and Gornitz, 1984; Delworth and Manabe, 1989; 
Oglesby and Erickson, 1989; Dirmeyer, 1994; Lau and Bua, 1998; Xue et al., 
2001, 2004; to name only a small fraction). Generally missing from these 
studies, however, is an analysis of the degree to which the experimental re- 
sults are model-dependent. Such model dependence can bias results tremen- 
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dously. Consider two hypothetical AGCMs, one in which the atmosphere 
responds strongly to anomalies in surface fluxes, which in turn respond to 
anomalies in land surface state, and one in which the atmosphere has an 
internal variability (chaotic dynamics) that overwhelms any signal from the 
land surface. (Note that in this paper, the term ’land surface’ refers to the 
combination of the vegetation canopy, the soil-atmosphere interface, and the 
top few meters or so of the soil, as typically modeled by AGCMs.) Experi- 
ments with these two AGCMs would lead to contradictory conclusions about 
the importance of properly initializing soil moisture in forecast simulations, 
about the degree to which deforestation affects climate, and perhaps even 
about the need for a realistic treatment of land surface processes in climate 
simulations. Contradictory results regarding land-atmosphere interaction do 
pervade the literature; see, for example, the broad range of results on defor- 
estation outlined by Hahmann and Dickinson (1997). 

The degree to which the atmosphere responds to anomalies in land sur- 
face state in a consistent manner, particularly at daily to seasonal timescales, 
is hereafter loosely referred to as the “land-atmosphere coupling strength”. 
Coupling strength cannot be determined a priori from a look at the model’s 
computer code. It is not explicitly prescribed or parameterized; it is rather a 
net result of complex interactions between numerous complex process param- 
eterizations in the AGCM, such as those for evapotranspiration, boundary 
layer development, and moist convection. Arguably, a shortcoming in the 
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analysis of the model-generated climate system is that this coupling strength, 
though a fundamental element of the system, is rarely examined closely and 
is almost never objectively quantified. The great majority of AGCM land- 
atmosphere interaction studies do not address the realism of the coupling 
strength implicit in the model used or how it compares with that in other 
models. 

An objective quantification and documentation of the coupling strength 
across a broad range of models would be valuable, if only to serve as a 
frame of reference when interpreting the experimental results of any par- 
ticular model. This objective documentation and quantification is indeed 
the goal of GLACE (for Global Land- Atmosphere Coupling Experiment). 
GLACE aims to show the extent to which coupling strength varies between 
models, and, more importantly, to characterize individual models as having 
a relatively strong, intermediate, or weak coupling, for later use in inter- 
preting various results obtained with those models. The range of coupling 
strengths uncovered by GLACE serves to quantify the uncertainty inherent 
in our understanding of land-atmosphere coupling and our ability to model 
it. 


1.2 Relationship to Pilot Study 

GLACE is a broad follow-on to the four-model intercomparison study of 
Koster et al. (2002), hereafter referred to as K02. K02 describes a numeri- 
cal experiment performed by four independent AGCM modeling groups, an 
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experiment that quantified, for each of the models, the degree to which pre- 
cipitation responds consistently to a prescribed, model-consistent time-series 
of land surface prognostic states. The chief result of K02 was a marked 
disparity in the coupling strengths of the four models. 

GLACE extends the K02 study substantially: 

- Participation From a Wider Range of Models. The intriguing inter- 

model variations discovered by K02 are presumably indicative of a 
broad range of coupling strengths implicit in today’s AGCMs. The 
goal of GLACE is to establish this range more precisely and (more 
importantly) to generate a comprehensive “table” of AGCM cou- 
pling strengths, a table that can help in the interpretation of the 
published results of a wide variety of climate models. 

- Separation of the Effects of “Fast” and “Slow” Reservoirs. The ex- 

perimental set-up used in K02 was limited; the prescribed diurnal 
surface temperature variations appeared to have had as much an 
effect on coupling strength as anything else. Since the initialization 
of surface temperature and water amounts in “fast” moisture reser- 
voirs (e.g., canopy interception) have little potential for prediction, 
particularly at subseasonal timescales and longer, the differences 
uncovered by K02 may have limited practical application. Of much 


5 



greater relevance to many land impact questions is whether some 
of the “slower” state variables - those variables with significant 
“memory” (in particular, soil moisture in the root zone and deeper 
reservoirs) - have an impact on the evolution of weather. This as- 
pect of coupling strength is a major focus of GLACE. 

- Effect on Air Temperature. K02 focused on how the land surface bound- 

ary affects the generation of precipitation. Also of interest is the 
control of the land boundary on air temperature fluctuations, par- 
ticularly when only root zone (and lower) soil moisture is prescribed. 
GLACE provides the means to address this issue. 

- Correction of Miscellaneous Technical Issues. Numerous technical prob- 

lems were encountered during the K02 study. Appropriate correc- 
tions are incorporated into the design of GLACE. The resulting 
model intercomparison is, as a result, cleaner. 

GLACE can indeed claim participation from a wider range of models. 
The experiment was offered to the community in early 2003. Over the course 
of that year, twelve AGCM groups performed the experiment and submitted 
results for processing. 
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1.3 Focus of Paper 


The purpose of the present paper is twofold. First, it thoroughly describes 
and contrasts the inherent coupling strengths of the twelve participating 
models. It thus provides a “snapshot” of the current state of land-atmosphere 
modeling, with emphasis (unlike K02) on the impacts of the slow reservoirs 
relevant to seasonal prediction. Second, and perhaps more important, it 
provides a full set of instructions for performing the experiments. This will 
allow additional models or future versions of the participating models to 
repeat them at will and immediately place their model’s behavior in the 
context of the behaviors documented here. A companion paper (Guo et al., 
this issue) examines the model-to-model differences in coupling strength - 
and the spatial variations in coupling strength seen within a given model - 
in the context of parameterization differences and differences in climatological 
and hydrological regime. 

Neither paper, however, addresses the realism of simulated coupling strength, 
primarily because direct measurements of land-atmosphere interaction at 
large scales do not exist. The identification of the proper measurements 
to make and their subsequent collection and analysis would clearly advance 
the study of this interaction. Potential local assessments of coupling strength 
and indirect large-scale evaluations are reserved for a future study. 

In the present paper, the experimental design is described in section 
2, with technical details relegated to an appendix. Section 3 provides an 
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overview of the participating models. Section 4 presents the basic results, 
and Section 5 provides a look at where on the globe the models tend to agree. 

2 Experimental Design 

2.1 Overall Framework 

The GLACE experiments consist of three separate 16-member ensembles of 
AGCM simulations, each simulation covering the period June 1 - August 
31. In CPU terms, this is equivalent to a single 12-year AGCM simulation. 
The overall design of the experiments is illustrated in Figure 1, and the run 
specifications are summarized in Table 1. 

The first ensemble, called Ensemble W, is essentially a standard set of 
AGCM simulations with prescribed sea surface temperatures. The only un- 
usual aspect of this ensemble is that in one of the simulations, chosen ran- 
domly but referred to here as “Wl” for convenience, all land prognostic vari- 
ables are recorded into a special data file at every time step. (See top panel of 
Figure 1.) The special data file is hereafter referred to as Wl-STATES. The 
recorded prognostic variables include soil moisture contents at all vertical lev- 
els, temperatures at all vertical levels, canopy interception reservoir content, 
and various variables characterizing snow, if snow is present. One global field 
is recorded per state variable per time step. K02, by the way, demonstrated 
that the choice of the ensemble member used to write into Wl-STATES is 
unimportant (at least for the one model examined); on a global average, any 
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ensemble member should produce approximately the same results in the later 
parts of the experiment. 

The second part of the experiment consists of another 16-member ensem- 
ble of 3-month simulations, using the same prescribed SSTs and the same 16 
sets of atmospheric initial conditions. In this ensemble (hereafter referred to 
as Ensemble R) , all member simulations are forced to maintain precisely the 
same time series of (geographically-varying) land surface states - namely, the 
states generated in simulation Wl. If for example, simulation W1 produced 
a very wet soil in southern France on July 27, then the atmosphere in every 
simulation of ensemble R is forced to feel the same very wet soil in southern 
France on July 27. This effect is achieved by discarding, at every time step 
of every R simulation, the updated values of all land surface prognostic vari- 
ables and then replacing them with the corresponding values for that time 
step from W1_STATES. (See middle panel of Figure 1.) 

The final part of the experiment, referred to as ensemble S, is equivalent 
to ensemble R except that only a small subset of the land surface prognostic 
variables are reset at each time step, as illustrated in the bottom panel of 
Figure 1. In particular, only soil moistures corresponding to soil layers with 
centers 5 cm or more below the surface are reset from W1.STATES. The 
other variables (e.g., temperatures, canopy interception contents, and soil 
moisture in a thin surface layer, if such a layer exists) are allowed to evolve 
freely, as they did in ensemble W. Most of the analysis in this paper and the 
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companion paper will focus on ensemble S, since it isolates and quantifies 
the impact of a relatively predictable state (deep soil moisture, a state with 
significant memory [Koster and Suarez, 2001]) on the evolution of weather. 

SST boundary conditions for all of the integrations correspond, as much 
as possible, to the period June-August 1994. This year was chosen because 
neither El Nino nor La Nina conditions during the year are strong. Different 
SST datasets are available, but for consistency, modeling groups were asked 
to use the AMIP-2 SST dataset (Gleckler, 1996) if at all possible. The 
K02 study, by the way, suggests that the impact of the chosen SST field on 
overall land-atmosphere coupling strength is small, though the choice of the 
year may perhaps have some bearing on specific geographical details. 

Specific details of the experimental design, including rules for initializa- 
tion of the different ensemble members, are provided in Appendix A. 

3 Participating Models 

Twelve AGCMs participated in the experiment. They are labeled as follows: 
BMRC, for the Bureau of Meteorology Research Centre in Melbourne, Aus- 
tralia; CCCma, for the Meteorological Service in Toronto, Canada; CCSR/NIES, 
for the Center for Climate System Research/National Institute for Environ- 
mental Studies in Tokyo, Japan; COLA, for the Center for Ocean-Land- 
Atmosphere Studies in Calverton, Maryland, United States; CSIRO-CC3, 
for the Commonwealth Scientific and Industrial Research Organization in 
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Aspendale, Victoria, Australia; GEOS-CRB, an AGCM used in the Climate 
and Radiation Branch at the Goddard Space Flight Center in Greenbelt, 
Maryland, United States; GFDL, for the Geophysical Fluid Dynamics Lab- 
oratory in Princeton, New Jersey, United States; HadAM3, an AGCM used 
at the Hadley Centre for Climate Prediction and Research in Exeter, United 
Kingdom; CAM3, for the Community Atmosphere Model (version 3) used at 
the National Center for Atmospheric Research in Boulder, Colorado, United 
States; GFS/OSU, which stands for “Global Forecast System model coupled 
to the Oregon State University land surface model”, used at the National 
Center for Environmental Prediction in Camp Springs, Maryland, United 
States; NSIPP, for the NASA Seasonal-to-Interannual Prediction Project, 
now part of the Global Modeling and Assimilation Office at the Goddard 
Space Flight Center in Greenbelt, Maryland, United States; and UCLA, for 
the University of California at Los Angeles, United States. Table 2 lists 
important details regarding the implementation of each of these models. 

4 Results 

4.1 Precipitation 

Using the Q diagnostic defined by K02, we examine here the land surface’s 
control on “synoptic-scale” precipitation variability, that is, the variability of 
precipitation on timescales of about a week. First, we aggregate the precip- 
itation output from each simulation into time-series of 6-day totals. Given 
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that the simulations are 92 days long and that we ignore the first 8 days to 
avoid problems associated with initial “shocks” to the modeled atmosphere, 
each simulation provides a time series, P(t), consisting of fourteen 6-day to- 
tals. For a given ensemble, which consists of 16 simulations, the standard 
deviation of precipitation, crp, is computed across the resulting 224 6-day 
totals. (The choice of 6 days for the time-aggregation is arbitrary; other 
choices give similar results.) 

Next, we compute the ensemble mean time series, P(t ): 

P(t) = iug.PiM, (i) 


where i represents the index of the ensemble member. The fourteen values 
in P(t) are used to compute the standard deviation of the ensemble mean 
time series, ap. Finally, ap and ap are combined into the diagnostic f Ip: 


p = 


16 a 2 p — o 2 P 
15cTp 


( 2 ) 


Qp measures the degree to which the sixteen precipitation time series gen- 
erated by the ensemble members are similar. In essence, f2p is equivalent to 


the ratio of explained precipitation variance to total precipitation variance; 
it varies from (approximately) 0 to 1, with higher values implying a greater 
contribution of boundary and initial conditions (and thus a lesser contribu- 
tion of atmospheric chaos) to the evolution of precipitation in a given AGCM. 
Notice that if all simulations produced precisely the same precipitation time 
series, implying no chaotic contribution, ap would be identical to crp, and Vt P 
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would be exactly one. K02 provides a graphic interpretation of the meaning 
of the Op diagnostic (see Figure 2 of K02). 

Figure 2 of the present paper shows the global fields of 0 p(W) (that 
is, O from the W ensemble) for all 12 AGCM-LSM combinations. Land 
states are not prescribed in ensemble W; thus, Op (IF) reflects the extent to 
which low frequency seasonal variations, as induced by the time variations 
of imposed boundary conditions and forcing alone (e.g., SST, vegetation 
structure, and solar declination) , lead to strong coherence in the precipitation 
rates. (Note that while this coherence may be strengthened in ensemble 
W through land-atmosphere feedback, the ultimate source of the coherence 
lies in the prescribed boundary conditions and forcing.) The high values of 
Qp(W) tend to be clustered in the tropics (where the ITCZ is migrating) 
and in a few midlatitude regions, such as eastern and southern Europe. An 
example of a model’s behavior at a grid cell with high Q, P (W) is shown in 
Figure 3. Plotted in the figure are sixteen time series of precipitation, one 
for each of the ensemble W simulations produced by CCSR/NIES over a 
grid cell in equatorial Africa. The same strong seasonality pervades each 
ensemble member, leading to a high synoptic-scale precipitation coherence 
over the duration of the simulation and thus to a high value (0.59) of Q, P (W). 

To quantify land-atmosphere coupling strength, we note that in ensem- 
ble R, the explained variance - the coherence of precipitation between the 
ensemble members - has two distinct sources: (a) the prescribed land vari- 
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ables and (b) the background seasonal behavior that contributes to £lp(W), 
as exemplified in Figure 3. Thus, subtracting fi, P for ensemble W from that 
for ensemble R should isolate the impact of prescribed land variables on 
the synoptic-scale precipitation variance. We use the difference in coherence 
f lp(R) — Qp(W) to measure land-atmosphere coupling strength associated 
with the prescription of all land variables. At a single grid cell, an differ- 
ence of 0.06 is significant at the 95% confidence level. 

Of course, if fJp(W) is already close to 1, the impact of land conditions 
will necessarily be small. This is not a major issue, however; the maximum 
of fip(W) across the different models over non-ice land points lies below 
0.8, and f2p(W) generally falls far below this maximum. For some models 
(GFS, BMRC), fip(W) values are small across the globe, implying that their 
“background” low frequency precipitation variations are quite small, at least 
for the 92-day period considered. 

Figure 4 shows the global fields of £lp(R) — £lp(W) for all 12 AGCM- 
LSM combinations. This figure is, in effect, a more comprehensive version 
of Figure 3 in K02. As in that earlier figure, Figure 4 shows a wide in- 
termodel disparity in the degree to which the atmosphere responds to the 
imposed land surface anomalies. Some models have relatively high values 
off l P (R) - ftp(W) (e.g., GFDL, NSIPP, CAM3, COLA, CSIRO), and oth- 
ers show relatively low values (e.g., HadAM3, BMRC, GFS, GEOS-CRB). 
Generally, however, f 'lp(R) — flp(W) is small in southern hemisphere mid- 
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latitudes and in deserts, presumably because the low mean evaporation rates 
imply little variability in the surface energy balance. The low evaporation 
rates in the summer hemisphere presumably reflect wintertime conditions; a 
repeat of the experiments for austral summer could prove useful. 

While the patterns and magnitudes shown in Figure 4 are intriguing, they 
may be of largely academic interest, since they may be controlled mostly by 
“fast” land surface prognostic variables - variables having little temporal 
memory and which cannot be used for prediction. In contrast, ensemble S 
focuses on the “slow” subsurface moisture variables that can contribute to 
prediction. Figure 5 shows the global fields of Qp(S) — Q,p(W) for all 12 
AGCM-LSM combinations. In analogy to Qp(R) — fip(W), Qp(S') — flp(W) 
isolates the contribution of prescribed subsurface soil moisture to precipita- 
tion variability - to the evolution of precipitation on synoptic timescales. 

As in the comparisons of Q,p(R) — fip(W), a strong diversity of coupling 
strength is seen amongst the AGCMs. Some models (e.g., GFDL, CAM3, 
NSIPP, CCCma) have a distinct “blocky” structure associated with their 
flp (5) — fip(W) values - large patches of relatively high Clp(S) — Clp(W) 
can be seen, for example, in central North America and the Sahel. Other 
models (e.g., CCSR/NIES, HadAM3, BMRC, GFS) have relatively few such 
structures; for the most part, small values of Qp(S) — flp(W) are scattered 
randomly across the globe. Note that a certain amount of agreement is seen 
in the positioning of the £lp(S) — flp(W) structures that do appear. This 
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will be discussed further in section 5. 

Most models even show some negative values of Q P (S) — ftp(VF). The 
reasons are unclear, but the highly infrequent negative values may have oc- 
curred by chance - according to monte carlo analysis, under an assumption 
of independent rainfall amounts in consecutive 6-day periods, a difference 
of either -0.1 or 0.1 is statistically significant at the 99.6% level, so a false 
negative may occur in about 0.4% of the grid cells plotted. Differences of 
-0.05 and 0.05 are statistically significant at the 92% level. (Relaxing the 
assumption of independence makes the occurrence of spurious negative val- 
ues slightly more likely.) Note that for the £l P (R) — Q P (W) field, negative 
values may have a different source; the specification of land states may have 
led to artifically large vertical gradients between the land surface and the 
free-running atmosphere, causing instabilities and unrealistic fluxes in the 
integrations and thus abnormal model behavior (Reale et al. 2002). 

For ease of comparison, Figure 6 shows the values of Op(lT), Q P (R) — 
Q P (W) and Q, P (S) — Q P (W) for each model, averaged across all non-ice 
land points. Note the different scale for the bottom plot; the specification 
of subsurface soil moisture has a much smaller impact on precipitation’s 
synoptic-scale variability than does either the background seasonality (top 
plot) or the specification of “fast” variables (middle plot). Though the num- 
bers for Qp(S’) — Qp(fF) appear small, we must keep in mind that the global 
averaging will hide any reasonably large Q P (S) — fip(VF) values that appear 
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regionally. In certain regions, subsurface soil moisture can have a significant 
impact on rainfall - and can thus be useful for seasonal prediction - even if 
the globally averaged impact appears small. 

The model diversity seen in the histograms reflects that seen in the maps. 
On a global average basis, the coupling strength associated with all land reser- 
voirs (middle plot of Figure 6, showing Clp(R) — Cl P (W)) is more than four 
times higher in some models (e.g., CAM3, NSIPP) than it is in some other 
models (e.g., BMRC, HadAM3, GEOS-CRB). The impact of subsurface soil 
moisture on the evolution of precipitation (bottom plot) is just as model- 
dependent; the more strongly coupled models (GFDL, CAM3, NSIPP, CC- 
Cma) stand out distinctly from the more weakly coupled models (HadAM3, 
BMRC, GFS). 

4.2 Surface Air Temperature 

The global fields of the Cl difference for surface air temperature are pre- 
sented in Figure 7. The temperature analysis focuses on ensemble S, since 
temperatures in ensemble R are overly influenced by the specified ground 
temperatures. In analogy to the precipitation analysis, the global fields of 
Clp(S) — fir(^) indicate a strong disparity in the control of subsurface soil 
moisture on the synoptic-scale variability of air temperature, with some mod- 
els (e.g., GFDL, HadAM3, CCCma, CSIRO) showing a high degree of control 
and others (e.g., COLA, BMRC, GFS, CCSR/NIES) showing a much more 
limited control. Also in analogy to the precipitation analysis, some regions 
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of coupling (e.g., the Sahel, northeastern China, and south-central North 
America) tend to show up in many of the models. 

The top panel of Figure 8 shows the global average (over non-ice land 
points) of Ot(VF). The averages are much larger than those for Qp(VF), 
presumably because of the strong background seasonal temperature cycle 
within each model. The middle panel of Figure 8 shows the average of 
Qt(S) — £Lt(W) across non-ice land squares for each model. Comparing 
this figure with the bottom panel of Figure 6 reveals two important things. 
First and foremost, the land’s control on air temperature is generally much 
larger than its control over precipitation. This is not surprising given that 
evaporation, through latent cooling and the associated impact on sensible 
heat flux, has a stronger connection to near-surface temperature than to pre- 
cipitation, which is produced at higher atmospheric levels and thus depends 
in part on convection and boundary layer formulations. (This is addressed 
in more detail in the companion paper [Guo et al., this issue].) Second, a 
low control on precipitation relative to other models does not imply a rela- 
tively low control on air temperature as well. HadAM3, for example, has a 
reasonably high average CIt(S) — fir(VF) values despite its low average value 
of Q P {S)-n P (W). 
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5 “Hot Spots” of Land- Atmosphere Coupling 

Figure 5 does suggest some similarity between models in precipitation’s re- 
sponse to soil moisture, at least in terms of geographical distribution. Several 
models, for example, place relatively high values of ttp(S) — Qp(W) in the 
Sahel and in central North America. Some intermodel similarity is also seen 
in the Qt(S) — fir(W0 fields. The GLACE experiment has a noteworthy 
strength: it provides a unique chance to quantify multi-model “agreement” 
in the locations of land moisture impact on the atmosphere. It can provide a 
more robust estimate of where the coupling is relatively strong, an estimate 
that is less subject to the quirks or deficiencies of any one particular model. 

This strength of GLACE motivated a recent paper (Koster et al., 2004) 
highlighting these “hot spots” for precipitation, i.e., identifying the areas 
which, for many of the models, the land-atmosphere coupling strength is rel- 
atively large. Plotted in that paper was the global field of Qp(S) — Qp(lT), 
averaged across all participating models. A slightly different version of the 
plot is shown in the top panel of Figure 9; the version is different because 
here, to maintain consistency with the rest of our 2-part paper, statistics are 
computed on the precipitation values themselves rather than on their natu- 
ral logarithms. (Although performing statistics on the natural logarithms of 
precipitation is a common and useful practice in hydrology and meteorology, 
since it reduces noise associated with high rainfall amounts, it produces tech- 
nical problems for some of our analyses.) To produce the plot, the results 
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from each model were disaggregated to the same very fine grid, one with a 
resolution of 0.5° x 0.5°. Disaggregation was performed in the simplest way 
possible. Each 0.5° x 0.5° grid cell lies wholly or mostly within a coarse grid 
cell of a given model. The precipitation rate assigned to the fine grid cell 
was that which applies to the coarse grid cell containing it. 

The top panel of Figure 9 shows that hot spots appear in the central 
Great Plains of North America, northern India, the Sahel, equatorial Africa, 
and a few additional regions. Because the logarithms of rainfall amounts 
are not used, however, the magnitudes of the plotted coupling strengths are 
slightly reduced relative to those in Koster et al. (2004). 

Note that a strict arithmetical average across the twelve models was used 
to generate the figure. An alternative approach would be to give added weight 
to the models with more realistic climate. The bottom panel of the figure 
shows the results of one such weighted calculation. In the approach used 
here, the Qp(S) — f2p(W) values are averaged across only the eight models 
that, at a given grid cell, best reproduce the observed climatological average 
precipitation for June through August (from GPCP [Huffmann et al, 1997]). 
Thus, a different set of eight models may contribute to the plotted average 
at adjacent locations. This approach is limited in scope; indeed, all possible 
weighting approaches are necessarily imperfect. The two chief deficiencies 
of the weighting applied here are that (i) rainfall rates used to evaluate the 
“realism” of a model may reflect the year chosen for the SST boundary 
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condition, whereas the climatological average for observations represents a 
mean over many years, and (ii) a model may have a realistic mean climatology 
but poor variability characteristics, and vice-versa. Nevertheless, the results 
of the exercise are illuminating. While the positions of the hot spots are 
similar to those in the top panel, the magnitudes of the averages have, in 
general, increased. In other words, by focusing on the models that appear 
to be more realistic in terms of precipitation climatology, we have increased 
the derived average coupling strength. 

The equivalent two maps for air temperature are plotted in Figure 10. 
The top panel shows £It(S) — Qt(W) averaged over all the models, and 
the bottom panel shows the “weighted average” result, again an eight-model 
average based on the realism of simulated precipitation. The results from 
both maps suggest strong synoptic-scale coupling for temperature in the 
Sahel, the central Great Plains of North America, India, and (in contrast 
to precipitation) eastern Asia. Notice that the average coupling strength is 
significantly larger than that for precipitation. 

Again, direct estimates of coupling strength from observations do not ex- 
ist, and the coupling results from any one model may simply reflect the pe- 
culiarities of that model. The multi-model averaging procedure, though sub- 
ject to deficiencies shared by multiple models, and though unable to generate 
quantitative estimates of reliability, still provides what is probably the best 
estimate possible for land-atmosphere coupling strength in the real world. 


21 


6 Summary 

In nature, rainfall certainly affects soil moisture, and soil moisture may affect 
rainfall. As part of the GLACE project, a number of AGCM groups have 
performed a numerical experiment designed to isolate the latter direction of 
causality. Through GLACE, we quantify the impact of land conditions on the 
evolution of precipitation and temperature in boreal summer in each of the 
models, and we compare in detail the differences in this “coupling strength” 
between the models. 

This paper has two main functions: (i) it describes the GLACE experi- 
ment with enough detail to allow its execution in the future by any modeling 
group, and (ii) it documents the range of coupling strengths implicit now 
in the participating models, so that any future group can put their results 
immediately into context. The range of coupling strengths uncovered by 
GLACE is indeed large, as indicated by Figures 4 through 8. We empha- 
size again that this intermodel disparity is not a trivial result, since coupling 
strength is not an explicitly defined quantity in the AGCMs - it is rather a 
complex product of many interacting model parameterizations. Most mod- 
elers have little notion of the degree of land-atmosphere coupling implicit in 
their models. The GLACE experiment provides, for the first time, an es- 
tablished methodology for its computation. Being able to compare a given 
model’s coupling strength to that of other models is critical for interpreting, 
for example, land use impact experiments or precipitation forecasts based on 
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soil water initialization. 

A side benefit of the GLACE experiments is the determination of multi- 
model “hot spots” of land-atmosphere coupling - regions that, according to 
several AGCMs, have a relatively high coupling strength. Figure 9 shows, 
for example, that the Sahel and the Great Plains of the U.S. are hot spots 
of coupling for precipitation at synoptic timescales. The multi-model nature 
of this result gives it added validity; either several models are wrong in a 
similar way, or these are indeed regions of strong coupling in the real world. 

Two questions naturally arise from this study. First, what causes the ge- 
ographical variations in coupling strength seen for a given model in Figures 
5 and 7? Second, what causes the model-to-model differences in coupling 
strength, as summarized by the histogram plots? The answers certainly re- 
late to differences in the parameterizations employed by the models and to 
differences in the simulated climates - some hydroclimatological regimes are 
presumably more amenable to coupling than others. Part 2 of this series of 
papers (Guo et al., this issue) addresses these two questions in detail. 


Appendix A: Details of Experimental design 


A.l Model-Specific Aspects of Experiment 
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The spatial resolution and the time step used necessarily varies amongst 
the participating AGCMs. Each group used a resolution typical for their 
model. Each group also applied their own strategy for writing out the prog- 
nostic variables in Ensemble W and for reading them in Ensembles R and S. 

A. 2 Initialization of Ensemble Members 

The members of an AGCM ensemble typically differ only in their initial 
atmospheric and land surface conditions. The approach for assigning the 
initial conditions is not strictly specified by GLACE; the only requirement 
is that the initial conditions be fully consistent with the AGCM being used. 
They are not allowed to be imported from some other model. 

Several approaches for initializing land and atmosphere states are possi- 
ble; they are listed in order of preference below. (That is, approach (a) is 
preferred most.) The key is to produce sets of initial conditions that sample 
the full range of possible land and atmosphere states. Initial land conditions 
between ensemble members, for example, should not be allowed to be artifi- 
cially similar, as can happen through the commonly used approach (e). 

(a) Some groups have available an archived series of 16 or more parallel 
multidecade AMIP-type simulations (i.e., simulations using SSTs 
prescribed from observations) from which to extract 16 different 
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sets of land and atmosphere states for June 1, 1994. These states 
can be used to initialize the W, R, and S ensembles. If daily data 
from the 16+ parallel AMIP-type simulations are archived, then in 
effect ensemble W is already almost finished; only one more simula- 
tion - the one that writes the time step information to Wl -STATES 
- needs to be performed for that ensemble. 

(b) If the number of archived multidecade AMIP-type simulations is less 

than 16 but greater than 1, they can still be used, as long as the 
years from which the June 1 land and atmosphere states are ex- 
tracted belong to the set of “quiescent” years (i.e., years with little 
El Nino or La Nina signal). For the purposes of this experiment, 
these years are 1951, 1952, 1959, 1960, 1961, 1963, 1977, 1979, 1980, 
1981, 1986, 1990, and 1994, years for which the Nino3 anomaly has 
an absolute value less than 0.5 for the three months preceding the 
initialization date. A group, for example, may have four archived 
parallel AMIP simulations. Extracting restart files for June 1 of 
1977, 1979, 1990, and 1994 from each of the 4 simulations would 
give a total of 16 sets of initial states for the experiment. 

(c) A more tractable approach for many groups is to access restart files 

(initial conditions) from a preexisting single 16+ year simulation. 
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In particular, if such a simulation exists in which SSTs do not vary 
from year to year (i.e., they are set to seasonally- varying climato- 
logical values), then the land and atmosphere states produced on 1 
June in each of 16 years of the simulation can be used to initialize 
the 16 ensemble members. 

(d) If the only 16+ year simulation available is an AMIP-type simulation 

(one with interannually- varying SSTs), then the June 1 conditions 
determined for the different years in this simulation can be used 
to initialize the June-August 1994 simulations. With this approach 
(as with approach (b)), the calendar years for the AMIP simulation 
are forced to lose their meaning. For example, suppose the restart 
files produced by an AMIP-type simulation covering 1979-1994 are 
available. The 1 June 1979 atmosphere and land states can be used 
to initialize one member of ensemble W (and of ensembles R and 
S), the 1 June 1980 states can be used to initialize another ensemble 
member, and so on. 

(e) A common approach to assigning initial conditions to the different 

members of an ensemble is to run the AGCM for, say, 16 June days 
and write out the atmosphere and land states at the beginning of 
each day. Each daily set of fields would then be used as initial 
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conditions. This type of approach, however, is highly undesirable 
for the present experiment, since the land surface states would not 
have time to vary much during the short simulation - the initial land 
conditions amongst the different members of ensemble W would not 
represent the broad range of states the model is capable of achieving. 

Note that given the design of the experiments, the initialization of all 
land states for ensemble R and the deeper soil moisture states for ensemble 
S is actually irrelevant. Note also that in all cases, the atmosphere may feel 
a “shock” at the beginning of the R and S simulations, since initially it will 
not be in equilibrium with the prescribed surface state. K02 examined the 
effect of this shock on Op and concluded that it was small. Nevertheless, 
the first 8 days or so of each 3-month simulation is excluded from the data 
analyses, to avoid its effects. 

A. 3 Energy and Water Balance Considerations 

The design of ensembles R and S necessarily precludes the maintenance 
of a strict energy and water budget below the land-atmosphere interface. 
Note, however, that energy and water in the atmosphere and across the 
interface are still perfectly conserved; conservation of energy and water is 
only “neglected” within the land reservoirs themselves. Since these special- 
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ized experiments focus solely on the atmospheric response to land conditions 
through the interface, the lack of conservation below the interface is deemed 
acceptable. 

A. 4 Redundancy of Simulations 

If the initial conditions used by simulation W1 (the simulation that wrote 
out its state variables into file W1 .STATES) are also used to initialize one of 
the members of ensemble R (say, simulation “Rl”), then by the construct of 
the experiment, the weather (and thus the precipitation) generated in simu- 
lations W1 and Rl should be identical. The same holds true if Wl’s initial 
conditions are used to initialize a member of ensemble S. Modeling groups 
can, if they wish, take advantage of this redundancy by using simulation W1 
as a member of both the R and S ensembles. In other words, in reality only 
15 new simulations need to be performed for both the R and S ensembles. 
(Note that truncation errors may, in fact, allow simulation Rl or SI to di- 
verge from simulation Wl. These truncation effects are irrelevant; the point 
is that simulation Wl can properly serve as a member of both the R and S 
ensembles.) 
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Figure Captions 


Fig. 1 Basic design of the experiment, as performed by all participating mod- 
els. 

Fig. 2 Global distributions of f)p(bF) for the models participating in GLACE. 

Fig. 3 Time series of rainfall (one line for each ensemble member) at a grid 
cell with a high fip(bF). 

Fig. 4 Global distributions of f 1 P (R) — £l P (W) for the models participating 
in GLACE. 

Fig. 5 Global distributions of £l P (S) — fip(bF) for the models participating 
in GLACE. 

Fig. 6 Mean of fl P (W) (top ),Q P (R)-Q P {W) (middle), and Q P (S)- Q P {W) 
(bottom) across non-ice land grid cells for the models participating in 
GLACE. 

Fig. 7 Global distributions of ftp (5) — flr(bF) for the models participating 
in GLACE. 

Fig. 8 Top: fFr(bF) averaged over non-ice land points. Bottom: flr(S) — 
J2cr(bF) averaged over non-ice land points.) 

Fig. 9 Top: Average of fip(S')— fip(bF) across all 12 models models. Bottom: 
average of fl P (S) — £l P (W) across the eight models that, at a given grid 
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cell, reproduce most closely the observed mean JJA precipitation. 


Fig. 10 Top: Average of 0, T (S) — £It(W) across all twelve models. Bottom: 
average of £It{S) — Qt(W) across the eight models that, at a given grid 
cell, reproduce most closely the observed mean JJA precipitation. 
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Ensemble 

Identifier 

# of Simulations 
in Ensemble 

Period Covered by 
Each Simulation 

Description 

W 

16 

June 1 - August 31, 1994 

Standard AGCM 
simulations with 
fully interactive 
land surface model. 

R 

16 

June 1 - August 31, 1994 

As W, except 
all land state 
variables replaced 
at every time step, 
from values in file. 

S 

16 

June 1 - August 31, 1994 

As W, except 
root zone and 
lower soil moisture 
variables replaced 
at every time step, 
from values in file. 


Table 1: Brief summary of GLACE ensembles. 




Model 


Resolution 


Prognostic variables 
set for R ensemble 


BMRC (Zhong et 
al. 2001; Colman 
etal. 2001; Des- 
borough, 1999; 
Desborough et al. 
2001 ) 


T47 


surface temperature; 
soil temperature; 
available moisture in 
root zone; canopy 
moisture storage; 
snow 


CCCma (McFar- 
lane et al. 1992; 
Boeretal. 1992; 
Verseghy 1991, 
2000; Verseghy et 
al, 1993) 


T32: 
3.75° x 
3.75° 


canopy temperature 
and water, snow 
temperature, depth, 
density, and albedo, 
soil temperature and 
moisture 


CCSR 
(Numaguti 1993; 
Numaguti et al. 
1997; Nozawa et 
al. 2001) 


T42 


surface temperature; 
soil temperature and 
moisture; frozen soil 
moisture; canopy 
water and temper- 
ature; snow 


COLA (Kinteret 
al. 1997; Xue et 
al. 1991, 1996; 
Dirmeyer and 
Zeng 1999) 


T63 

1.875° 


soil temperature and 
moisture for 3 layers; 
canopy interception; 
snow 


CSIRO-CC3 
(McGregor and 
Dix 2001; 
McGregor 1996; 
Kowalczyk et al. 
1994) 


2° x 2° 


soil temperature, 
moisture, ice; snow 
variables; canopy 
water reservoir 


GEOS (Conaty 
et al. 2001; Sud 
and Walker 
1999a,b; Mocko 
and Sud, 2001) 


2.5° x 2° 


soil moisture; ground 
temperature; canopy 
temperature and 
water amount; snow 
temperature, amount 
and density 


Prognostic 
variables set 
forS 
ensemble 

Air temperature 
variable 

Initializa- 
tion method 

available 
moisture in 
root zone 

air temperature at 
the lowest model 
layer (-50 m AGL) 

typeC 

Soil mois - 
ture,ice for 
all 3 layers 
(depth of top 
layer is 10 
cm) 

diagnosed 2m air 
temperature 

typeC 

soil moisture 
for all layers 
except the 
surface layer 

air temperature at 
the lowest model 
layer 

type D 

root zone 
and recharge 
layer mois - 
tures 

canopy air space 
temperature 

mix of type 
A and E. 

soil moisture 
at the 2 nd to 
6 th levels 

1 .8m screen air 
temperature 

hybrid of 
types D and 
E. 

soil moisture 
at root zone 
and deep 
recharge 
layers 

interpolated IO- 
meter above the 
ground 

modified 
version of 
type D 


Table 2: AGCM-LSM combinations participating in GLACE. 
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Model 

Resolution 

Prognostic variables 
set for R ensemble 

Prognostic 
variables set 
forS 
ensemble 

Air temperature 
variable 

Initializa- 
tion method 

GFDL (Milly and 
Shmakin 2002; 
GAMDT, 2004; 
but with different 
parameterization 
for boundary 
layer turbulence, 
prognostic clouds, 
and cumulus 
processes ) 


single column soil 
water; snow; soil 
temperature at 5 
vertical levels 

single 
column soil 
water 

interpolated 2m air 
temperature 

hybrid of 
types A and 
D. 

HadAM3 (Pope 
et al., 2000; Cox 
et al. 1999; Es- 
sery et al. 2003) 

3.75° x 
2.5° 

soil moisture and 
temperature; canopy 
water and snow at 
each tile 

soil moisture 
at all layers: 
top layer is 
lOcmdeep 

interpolated 1 .5m 
air temperature 

typeC 

CAM3 (Collins 
et al. 2004; Bonan 
et al. 2002; Ole- 
son et al. 2004) 

T42, 

(~2.8°x 

2.8°) 

snow properties; soil 
liquid water and ice 
contents; tempera- 
tures for soil, vegeta- 
tion, ground and lake; 
canopy water 

soil liquid 
water and 
ice contents 
below 5 cm 
depth 

air temperature at 2 
meters above the 
apparent sink for 
sensible heat 

typeC 

H 

T62, 

1.875° 

soil moisture and 
temperature at two 
layers; canopy water 
content; snow 

soil moisture 
in the 

second layer 

interpolated 2m air 
temperature 

type B 

NSIPP (Bacmeis- 
ter et al. 2000; 
Koster and 
Suarez, 1996) 

■ 

soil moisture, temper- 
ature, canopy inter- 
ception, and snow at 
each subgrid tile 

root zone 1 

and recharge 
layer 
moistures 

diagnosed 2m air 
temperature 

type A 

UCLA (Xue et 
al. 2001,2004) 

T42, 

2.5° x 2° 

soil moisture, tem- 
peratures; canopy air 
temperature, 
interception; snow. 

root zone 
and recharge 
layer soil 
moistures 

canopy air space 
temperature 

typeB 


Table 2: (cont.) 
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Part 1: Ensemble W (16 members) 


time step n+1 


time step n 


Step forward the 
coupled AGCM-LSM 


\(for W1) 


Write the values 
of the land surface 
prognostic variables 
Into file WI^STATES 


Step forward the 
coupled AGCM-LSM 


\(lpf W1) 


Write the values 
of the land surface 
prognostic variables 
into (He W1_STATES 


Part 2: Ensemble R (16 members) 



Part 3: Ensemble S 



Figure 1: Basic design of the experiment, as performed by all participating 
models. 
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Figure 2: Global distributions of Qp(V 
GLACE. 


44 



for the models participating in 




Precipitation: CCSR/NIES (lat = 1.395 degrees S; Ion = 25.312 E) 

' ' 1 ' ' 1 1 — ' , 1 t r 



Figure 3: Time series of rainfall (one line for each ensemble member) at a 
grid cell with a high fip(lF). 
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fl P (R)-fip(W) 



Figure 4: Global distributions of Qp(R)—S}p(W) for the models participating 
in GLACE. 
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fi P (s)-n„(w) 



Figure 5: Global distributions of Q P (S) -Q P (W) for the models participating 
in GLACE. 






Figure 6: Mean of fip(l'F) (top), £lp(R) — Q,p(W) (middle), and Qp(S) — 
Qp(fF) (bottom) across non-ice land grid cells for the models participating 
in GLACE. 
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Figure 7: Global distributions of £It(S) — f2r(W) for the models participating 
in GLACE. 



(a) Average of ft T (W): 




3 ) 

CD 


Figure 8: Top: iFr(W) averaged over non-ice land points. Bottom: &t(S) 
0.t{W) averaged over non-ice land points. 
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Figure 9: Top: Average of Q P (S) - Q P (W) across all 12 models models. 
Bottom: average of Q P (S) - fi P (W) across the eight models that, at a given 
grid cell, reproduce most closely the observed mean JJA precipitation. 
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Figure 10: Top: Average of Qt(S)— Q.t(W) across all twelve models. Bottom: 
average of £It(S) — Q.t(W) across the eight models that, at a given grid cell, 
reproduce most closely the observed mean JJA precipitation. 
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Popular Summary: 


“GLACE: The Global Land- Atmosphere Coupling Experiment. Part 1: Overview, 
and Part 2: Analysis”, by R. Koster, P. Dirmeyer, Z. Guo, and 21 others. (Submitted 
to Journal of Hydrometeorology) 


A wetter-than-usual soil may lead to higher-than usual evaporation, which in turn 
may lead to increased precipitation. This soil moisture - precipitation connection, if 
verified and properly utilized, would contribute significantly to seasonal forecasting 
efforts. Seasonal forecasters could then take advantage of the fact that initialized soil 
moisture anomalies can persist for months. 

The problem with verifying the soil moisture - precipitation connection with ob- 
servational data is that the required data on the large scale do not exist and are 
logistically impossible to obtain. Climatologists have thus relied instead on modeling 
studies to quantify the connection. These modeling studies have their own limitations, 
however; most notably, the results can be strongly model dependent. 

To forward our understanding of the soil moisture - precipitation connection, and 
in particular to address the question of model dependence in published results, the 
authors have operated GLACE, an international intercomparison project designed to 
quantify the strength of the soil moisture - precipitation connection (the “coupling 
strength” ) across a broad range of atmospheric general circulation models. Through 
GLACE, we find that the different models do indeed show a broad disparity in cou- 
pling strength distribution. GLACE, however, also provides an intriguing result. 
Despite the intermodel disparity, certain areas of the Earth show a large coupling 
strength in many models, suggesting that the existence of signifcant coupling strength 
in these areas is not so model-dependent. Given the lack of observational data, such 
a multi-model determination of areas with strong coupling strength is arguably the 
best estimate of such areas attainable by any method. 

Part 1 of the paper has two key goals: to document the intermodel variability in 
coupling strength existing in models today, and to provide a full set of instructions 
for repeating the experiment, so that other groups can test their models and compare 
their results directly with those documented in the paper. Part 2 (with Z. Guo as 
lead author) delves into the “whys” of the intermodel disparity, explaining in general 
terms what controls the coupling strengths of the different models. Together, the 
papers document the key results of a scientifically productive project. 



